Cosmological and astrophysical limits on 
brane fluctuations 



J. A. R. Cembranos, A. Dobado and A. L. Maroto 

Departamento de Fisica Teorica, 
Universidad Complutense de Madrid, 28040 Madrid, Spain 

ABSTRACT 

We consider a general brane-world model parametrized by the brane tension 
scale / and the branon mass M. For low tension compared to the funda- 
mental gravitational scale, we calculate the relic branon abundance and its 
contribution to the cosmological dark matter. We compare this result with 
the current observational limits on the total and hot dark matter energy 
densities and derive the corresponding bounds on / and M. Using the nu- 
cleosynthesis bounds on the number of relativistic species, we also set a limit 
on the number of light branons in terms of the brane tension. Finally, we es- 
timate the bounds coming from the energy loss rate in supernovae explosions 
due to massive branon emission. 



PACS: 95.35. +d, 11.25-w, ll.lOKk 



1 Introduction 



The increasing observational precision is making cosmology a useful tool in 
probing certain properties of particle physics theories beyond the Standard 
Model (SM). The limits on the neutrino masses and the number of neutrino 
families are probably the most clear examples [1, 2]. In some cases, the 
cosmological bounds are complementary to those obtained from colliders ex- 
periments and therefore, the combination of both allows us to restrict the 
parameter space of a theory in a more efficient way. 

There are two main ways in which cosmology can help. On one hand 
we have the relative abundances of the light elements, which is one of the 
most precise predictions of the standard cosmological model. Indeed, the 
calculations are very sensitive to certain cosmological parameters and thus 
for instance, the production of 4 He increases with the rate of expansion of the 
universe H . A succesful nucleosynthesis requires that H should not deviate 
from its standard value more than around a 10% during that epoch. Since 
H depends on the effective number of relativistic degrees of freedom g e ff(T) 
at a given temperature T, the above constraint translates into a bound on 
9eff(T nuc ), where T nuc ~ lMeV. Apart from the number of particle species, 
cosmology also sets a limit on their energy density. For particles which are 
non relativistic at present, there exists an upper bound given by the measured 
dark matter density Vt M = 0.23 ± 0.08 at the 95% C.L [1]. 

On the other hand, stars also provide useful information on theories con- 
taining light and weakly interacting particles. The extreme opacity of ordi- 
nary matter to photons makes it very long the time it takes for a photon 
produced in the center of the star to reach its surface. Indeed, this fact 
explains the longevity of stars. However, particles like neutrinos or axions, 
still can be produced abundantly in nuclear reactions in the core of the star, 
but since they are weakly interacting, the rate at which they can carry away 
energy can be much larger. This was particularly evident in the 1987A su- 
pernova explosion, in which most of the energy was released in the form 
of neutrinos. Again, this fact can be used to set limits on the mass and 
couplings of the new particles. 

In this paper we will study the constraints that cosmology and astro- 
physics impose on the so called brane- world scenario (BWS), which is be- 
coming one of the most popular extensions of the SM. In these models, the 
Standard Model particles are bound to live on a three-dimensional brane em- 



1 



bedded in a higher dimensional (D = 4 + N) space-time, whereas gravity is 
able to propagate in the whole bulk space. The fundamental scale of gravity 
in D dimensions M D can be lower than the Planck scale M P . In the original 
proposal in [3], the main aim was to address the hierarchy problem, and for 
that reason the value of Mp was taken around the electroweak scale. How- 
ever more recently, brane cosmology models have been proposed in which 
Mjj has to be much larger than the TeV [4, 5] . In this work we will consider 
a general BWS with arbitrary fundamental scale M D . 

The existence of extra dimensions is responsible for the appearence of new 
fields on the brane. On one hand, we have the tower of Kaluza-Klein modes of 
fields propagating in the bulk space, i.e. the gravitons. On the other, since 
the brane has a finite tension / 4 , its fluctuations will be parametrized by 
some ir a fields called branons. These fields, in the case in which traslational 
invariance in the bulk space is an exact symmetry, can be understood as 
the massless Goldstone bosons arising from the spontaneous breaking of that 
symmetry induced by the presence of the brane [6, 7]. However, in the most 
general case, translational invariance will be explicitly broken and therefore 
we expect branons to be massive fields. 

It has been shown [8] that when branons are properly taken into account, 
the coupling of the SM particles to any bulk field is exponentially suppressed 
by a factor exp(-M|-^M| ) /(87r 2 / 4 )), where M KK is the mass of the corre- 
sponding KK mode. As a consequence, if the tension scale / is much smaller 
than the fundamental scale M D) i.e. / <C M D) the KK modes decouple from 
the SM particles. Therefore, for flexible enough branes, the only relevant 
degrees of freedom at low energies in the BWS are the SM particles and 
branons. 

The phenomenological implications of KK gravitons for colliders physics, 
cosmology and astrophysics have been studied in a series of papers (see [9] 
and references therein), and the corresponding limits on M D and/or the 
number of extra dimensions N have been obtained. In the case of branons, 
the potential signatures in colliders have been studied in the massless case in 
[10] and in the massive one in [11]. Limits from supernovae and modifications 
of Newton's law at small distances in the massless case were obtained in 
[12]. Moreover in [13] the interesting possibility that massive branons could 
account for the observed dark matter of the universe was studied in detail. 
The main aim of the paper is to analyse the cosmological and astrophysical 
limits on the BWS through the effects due to massive branons. To that end 
we will be assuming that the evolution of the universe is standard up to a 
temperature around /. Indeed, this is the case of realistic brane cosmology 
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models in five dimensions [4]. Also more recently, six-dimensional models 
have been proposed in which the Friedmann equation has the standard form 
for arbitrary temperature [5]. 

The paper is organized as follows: in section 2 we give a brief introduc- 
tion to the dynamics of massive branons. Section 3 contains a summary of 
the main steps used in the standard calculations of relic abundances gen- 
erated by the freeze-out phenomenon in an expanding universe. In section 
4, we give our results for the thermal averages of branon annihilation cross 
sections into SM particles. Section 5 is devoted to the study of the limits 
on the branon mass and on the brane tension scale from the cosmological 
dark matter abundances. In section 6, after reviewing the limits imposed by 
nucleosynthesis on the number of relativistic species, we apply them to the 
case of branons. Section 7 contains the calculation of the rate of energy loss 
from a supernova core in the form of branons and the corresponding limits 
on / and M. Finally, section 8 includes the main conclusions of the paper. 
In an Appendix we have included the explicit formulas for the creation and 
annihilation cross sections of branons pairs. 

2 The branon field 

In this section we will briefly review the main properties of massive brane 
fluctuations (see [7, 14, 11] for a more detailed description). We will con- 
sider a single-brane model in large extra dimensions. Our four-dimensional 
space-time M4 is embedded in a D-dimensional bulk space which, for sim- 
plicity, we will assume to be of the form M D = M 4 x B. The B space is 
a given N-dimensional compact manifold, so that D = 4 + N. The brane 
lies along M4 and we neglect its contribution to the bulk gravitational field. 
The coordinates parametrizing the points in Mjj will be denoted by (x^, y m ), 
where the different indices run as fj, — 0, 1, 2, 3 and m— 1, 2, N. The bulk 
space M D is endowed with a metric tensor which we will denote by Gmn, 
with signature (+, — , — ). For simplicity, we will consider the following 

ansatz: 



The position of the brane in the bulk can be parametrized as Y M = (x^, Y m (x)), 
with M = 0, . . . ,3 + N and where we have chosen the bulk coordinates so 
that the first four are identified with the space-time brane coordinates x^ . 




(1) 
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We assume the brane to be created at a certain point in B, i.e. Y m (x) = Y™ 
which corresponds to its ground state. We will also assume that B is a homo- 
geneous space, so that brane fluctuations can be written in terms of properly- 
normalized coordinates in the extra space: 7r Q (:r) = f 2 Y a (x), a = 1, . . . , N. 
The induced metric on the brane in its ground state is simply given by the 
four-dimensional components of the bulk space metric, i.e. g^ v = = G^. 
However, when brane excitations are present, the induced metric is given by 



9,u = d,Y M d u Y N G MN (x, Y(x)) = g, u (x, Y(x)) - d,Y m d u Y n ~g' mn (Y(x)) (2) 



The contribution of branons to the induced metric is then obtained expanding 
(2) around the ground state [7, 14, 11]: 



Branons are the mass eigenstates of the brane fluctuations in the extra-space 
directions. The branon mass matrix M a p is determined by the metric prop- 
erties of the bulk space and, in the absence of a general model for the bulk 
dynamics, we will consider its elements as free parameters (for an explicit 
construction see [15]). Therefore, branons are massless only in highly sym- 
metric cases [7, 14, 11]. 

Since in the limit in which gravity decouples Mo — > oo, branon fields 
still survive [16], branon effects can be studied independently of gravity. The 
mechanism responsible for the creation of the brane is in principle unknown, 
and therefore we will assume that the brane dynamics can be described by 
a low-energy effective action derived from the Nambu-Goto action [7]. Also, 
branon couplings to the SM fields can be obtained from the SM action defined 
on a curved background given by the induced metric (2), and expanding in 
branon fields. Thus, the complete action, up to second order in n fields, 
contains the SM terms, the kinetic term for the branons and the interaction 
terms between the SM particles and the branons: 



where Tg^^g^J) is the conserved energy-momentum tensor of the Standard 
Model evaluated in the background metric. 



(3) 
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Tsm — ~ g^^-SM + 2 



5CsM 



(5) 



It is interesting to note that under a parity transformation on the brane, 
the branon field changes sign if the number of spatial dimensions of the brane 
is odd, whereas it remains unchanged for even dimensions. Accordingly, 
branons on a 3-brane are pseudoscalar particles. This implies that if we 
want to preserve parity on the brane, terms in the effective Lagrangian with 
an odd number of branons would be forbidden. 

The quadratic expression in (4) is valid for any internal B space, regard- 
less of the particular form of the metric g' mn . In fact the low-energy effective 
lagrangian is model independent and is parametrized only by the number 
of branon fields, their masses and the brane tension. The dependence on 
the geometry of the extra dimensions will appear at higher orders. These 
effective couplings thus provide the necessary tools to compute cross sections 
and expected rates of events involving branons in terms of / and the branons 
masses only. 

From the previous expression, we see that since branons interact by pairs 
with the SM particles, they are necessarily stable. In addition, their cou- 
plings are suppressed by the brane tension / 4 , which means that they could 
be weakly interacting, and finally, according to our previous discussion, in 
general, they are expected to be massive. As a consequence their freeze-out 
temperature can be relatively high, which implies that their relic abundances 
can be cosmologically important. 

3 Relic branon abundances 

In order to calculate the thermal relic branon abundance, we will use the 
standard techniques given in [17, 18] in two limiting cases, either relativistic 
(hot) or non-relativistic (cold) branons at decoupling. In this section we will 
review the basic steps of the calculation method. 

The evolution of the number density n a of branons 7i a , a = 1, . . . , N with 
N the number of different types of branons, interacting with SM particles in 
an expanding universe is given by the Boltzmann equation: 
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where 



a A = J2°(^ a ^X) (7) 

x 

is the total annihilation cross section of branons into SM particles X summed 
over final states. The —3Hn a term, with H the Hubble parameter, takes into 
account the dilution of the number density due to the universe expansion. 
These are the only terms which could change the number density of branons. 
In fact, since branons are stable they do not decay into other particles and 
since they interact always by pairs the conversions like 7r a X — > ir a Y do not 
change their number. Notice that we are considering only the lowest order 
Lagrangian and assuming that all the branons have the same mass. This 
implies that each branon species evolves independently. Therefore in the 
following we will drop the a index. 

The thermal average (<jav) of the total annihilation cross section times 
the relative velocity is given by: 

1 r d 3 pi d 3 p 2 w(s) 

M = < J W?^f f {El)f ( (8) 

where: 



. . ^ ^ sa A / AM 2 . , 

w(s) = E x E 2 o A v rel = — Wl — (9) 

The Mandelstam variable s can be written in terms of the components of 
the four momenta of the two branons pi and p 2 as s = (pi + p 2 ) 2 = 2(M 2 + 
E\E 2 — |pi||P2| cos#). Assuming vanishing chemical potential, the branon 
distribution functions are: 

m = ^ do) 

with a = for Maxwell-Boltzmann and a = — 1 for Bose-Einstein. In the case 
of non-relativistic relics T 3M, the Maxwell-Boltzmann distribution is a 
good approximation and we will use it for simplicity instead of Bose-Einstein. 
Finally, the equilibrium abundance is given by: 

d 3 p 

n, 



From (4), the thermal average will include, to leading order, annihilations 
into all the SM particle-antiparticle pairs. If the universe temperature is 



6 



above the QCD phase transition (T > T c ), we consider annihilations into 
quark- ant iquark and gluons pairs. If T < T c we include annihilations into 
light hadrons. For the sake of definiteness we will take a critical temperature 
T c ~ 170 MeV and a Higgs mass m H ~ 125 GeV, although the final results 
are not very sensitive to the concrete value of these parameters. 

In order to solve the Boltzmann equation we introduce the new variables: 
x = M/T and Y = n/s with s the universe entropy density. We will assume 
that the total entropy of the universe is conserved, i.e. S = a 3 s = const, 
where a is the scale factor of the universe and we will make use of the Fried- 
mann equation: 



H 2 = T p (12) 

3M 2 H V ; 



8tt 

where the energy density in a radiation dominated universe is given by: 

P = geff(T)^T* (13) 
In a similar way, the entropy density reads: 

s = h ff (T) 2 -^T* (14) 

where g e ff(T) and h e ff(T) denote the effective number of relativistic degrees 
of freedom contributing to the energy density and the entropy density respec- 
tively at temperature T (T being the temperature of the photon background). 
Notice that for T > MeV we have h e ff ~ g e ff. Using these expressions we 
get: 

where we have ignored the possible derivative terms dh e ff/dT. 

The qualitative behaviour of the solution of this equation goes as follows: 
if the annihiliation rate defined as Ta = n eq (<jAv) is larger than the expansion 
rate of the universe if at a given x, then Y(x) ~ Y eq (x), i.e., the branon 
abundance follows the equilibrium abundances. However, since decreases 
with the temperature, it eventually becomes similar to H at some point 
x = Xf. From that time on branons are decoupled from the rest of matter 
or radiation in the universe and its abundance remains frozen, i.e. Y(x) ~ 
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Y eq (xf) for x > Xf. For relativistic (hot) particles, the equilibrium abundance 
reads: 

y ^^WkJM' <*« 3) < 16) 

whereas for cold relics: 

= S (if * s/2 ^^ <*» 3 > < 17 > 

We see that for hot branons the equilibrium abundance is not very sensitive to 
the value of x. In the case of cold relics however, Y eq decreases exponentially 
with the temperature, which implies that the sooner the decoupling occurs 
the larger the relic abundance. 

Let us first consider the simple case of hot branons. Since its equilibrium 
abundance depends on ij only through h e ff(xf), the relic abundance is not 
very sensitive to the exact time of decoupling. In this case, in order to 
calculate the decoupling temperature Tf = M/xf, it is a good approximation 
to use the condition Ta = H . From the explicit expression of the Hubble 
parameter in a radiation dominated universe we have: 

H(T f ) = 1.67^(7))^ = F A (T f ) (18) 

which can be solved explicitly for Tf, expanding T^Ty) for Tf S> M/3. Once 
we know Xf, the relic abundance today (Y^ ~ Y{xf)) is given by (16). From 
this expression we can obtain the current number density of branons and the 
corresponding energy density which is given by: 

tt Br h 2 = 7.83-l(T 2 , - ^ (19) 
h eff (x f )eV 

The calculation of the decoupling temperature in the case of cold branons 
is more involved. The well-known result is given by: 

t>= jwwj (20) 

where c ~ 0.5 is obtained from the numerical solution of the Boltzmann 
equation. This equation can be solved iteratively. The corresponding energy 
fraction reads: 

Q Br h 2 = 8.77 • 10- n GeV- 2 4^ ( £ ^rxjA (21) 

9eff \n=0 n + 1 / 
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where we have expanded (<jav) in powers of x 1 as: 



oc 



(a A v) = °n X 



(22) 



n=0 



Notice that in general, oc l/(a A v), i.e. the weaker the cross section the 
larger the relic abundance. This is the expected result, since, as commented 
before the sooner the decoupling occurs, the larger the relic abundance, and 
decoupling occurs earlier as we decrease the cross section. Therefore the 
cosmological bounds work in the opposite way as compared to those coming 
from colliders. Thus, a bound such as Qsr < 0(1) translates into a lower 
limit for the cross sections and not into an upper limit as those obtained 
from non observation in colliders. 

In the following we apply the previous formalism to obtain the relic abun- 
dance of branons Qsrh 2 , both when they are relativistic and non-relativistic 
at decoupling. For that purpose we need to evaluate the thermal averages 
(ctav) for the annihilation of branons into photons, massive W ± and Z gauge 
bosons, three massless neutrinos, charged leptons, quarks and gluons (or light 
hadrons) and a real scalar Higgs field, in terms of the brane tension / and 
the branon mass M. 

4 Branon annihilation cross sections: ther- 
mal average 

We give the results for the different channels contributing to the thermal 
average of the annihilation cross section {(Jav) of branons into SM particles. 
The explicit production and annihiliation cross section can be found in the 
Appendix. For cold relics, we have expanded the expressions for each particle 
species in powers of 1/x as follows: 



In the hot branons case, we give the results for the different contributions 
to the decay rate T A = n eq (aAv), where we have considered the ultrarrel- 
ativistic limit for the branons, i.e. M — 0. For massive SM particles, the 
final expressions cannot be given in closed form. Therefore, in this section, 
in order to show the high-temperature behaviour, we only give the results for 
fermions, gauge bosons and scalars in the limit in which their masses vanish. 
Also in this case, we have used the Bose-Einstein form as the equilibrium 
distribution. 




(23) 
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4.1 Dirac fermions 



x > 3 (Cold) 



c ° = T^ M H(^ 2 -^)\/l-§ (24) 



1 / 777 ^ 

ci = T^^ M H(67M 2 - 31™^!-^ (25) 
1 M 2 



c 2 = ^ o „~ 9 ,„ a J9 "" — (17408M 6 + 13331M 4 m%. 



7680tt 2 / 8 M 2 - m 2 



TP 



- 46606M 2 m 4 + i8927mjWl-^| (26) 

Notice that this expansion is not valid near SM particles thresholds, i.e. for 
branon masses close to some SM particle mass. In addition, since the c 
coefficient is different from zero, annihilation will mainly take place through 
s-wave. 

x < 3 (Hot) 

For massless fermions, we obtain: 

Rir 9 T 9 

Tf rac = ° F77^ + W (27) 



297675C(3)/ 8 

4.2 Massive gauge field 
a; > 3 (Cold) 

M 2 yjl-!j£ (4M 4 -4M 2 m z 2 + 3m z 4 ) 
C ° = 6Tf^ (28) 



Cl = 



M 2 yjl-^ (364 M 6 - 584 M 4 m z 2 + 349 M 2 m z A - 93 m z 6 ) 
768 / 8 (M 2 - m z 2 ) 7T 2 



TVf 2 «/i — mz 

c 2 = i — 5 — (415756 M 8 — 755844 M 6 m z 2 

30720 / 8 (M 2 — m z 2 ) n 2 ^ 



+ 356541 M 4 m z 4 - 76294 AT m z b + 56781 m z 

Again this expansion is not valid near SM particles thresholds, and the lead- 
ing term corresponds to the s-wave. 
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x < 3 (Hot) 

In the limit where T ^> m z , one can obtain the following expression: 

4.3 Massless gauge field 
x > 3 (Cold) 

c = (30) 
ci = 

68 M 6 



c 2 



15/ 8 7T 2 



In this case, c = Ci = and the leading term corresponds to rf-wave annihi- 
lation. 



x < 3 (Hot) 



1 

r ^ 297675 C (3)/ S + °W < M » 



4.4 Complex scalar field 

x > 3 (Cold) 

M 2 (2M 2 + m, 2 )Vl-^ 
C ° = 32]V (32) 

M 2 (2 M 2 + m$ 2 ) y/l-^r (182 M 4 — 115 M 2 m^ 2 - 31 m^) 
° l = 384 / 8 (M 2 — m$ 2 ) 7T 2 



_/\/f2 A _ m*£ 

c 2 = i — (92164 M 8 — 123556 M 6 m$ 2 

5120 / 8 (M 2 -m$ 2 ) 2 7r 2 V 

+ 12269 M 4 m$ 4 + 9754 M 2 m$ 6 + 6309 m$ 8 ) 

We find the same problems near SM particles thresholds. The dominant 
contribution is the s-wave. 
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x < 3 (Hot) 

For massless scalar, one can obtain the following expression: 

1 fivr 9 T 9 

For real scalar fields, the above results should be divided by two. 

Notice that for conformal matter the leading contribution is the ci-wave. 
This explains why in the massless limit for fermions, the leading contribution 
is no longer the s-wave, but the ci-wave, whereas for massless scalars or taking 
the mz — > limit for massive gauge bosons, the s- and p-waves survive. 

Concerning the validity of the above results, in order to avoid the men- 
tioned problems of the Taylor expansion near SM thresholds, we have taken 
branon masses sufficiently separated from SM particles masses where the 
usual treatment is adequate [17, 18]. Such treatment is known to introduce 
errors of the order of 10% in the relic abundances. In addition, coannihila- 
tion effects are absent in this case since there are no slightly heavier particles 
which eventually could decay into the lightest branon. 

5 Cosmological bounds from the dark matter 
energy density 

For cold branons, once we know the c n coefficients for the total cross section, 
we can compute the freeze-out value Xf from (20) and the relic contribution 
to the energy density of the universe flsrh 2 from (21) in terms of / and M. 
Imposing the observational limit on the total dark matter energy density 
from WMAP: Vt Br h 2 < 0.129 - 0.095 at the 95% C.L. which corresponds to 
Vl M = 0.23 ± 0.08 and h = 0.79 — 0.65 [1], we obtain the exclusion plots in 
Figs. 1 and 6. 

Notice that in Fig. 6 we have plotted the Xf — 3 curve, which limits the 
range of validity of the cold relic approximation. Therefore, the excluded 
region is that between the two curves. It is also important to note that for 
those values of the parameters on the solid line, branons would constitute all 
the dark matter in the universe. 

For hot branons, we have computed numerically the total annihilation 
rate into SM particles r^. Using equation (18), we can find the freeze-out 
temperature Tf in terms of the brane tension scale /. Approximately, the re- 
lation between the logarithms of these quantities is linear: log 10 (//lGeV) ~ 
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/ (GeV) 




M (GeV) 

Figure 1: Exclusion plot from cold relics abundance. The thick lines corre- 
spond to Vlsrh 2 = 0.129 — 0.095 for N = 1 and N = 7. Therefore the areas 
above the curves corresponding to Qsrh 2 > 0.129 are excluded. The striped 
region contains the curves corresponding to 1 < iV < 7. 

(7/8) log 10 (T//l GeV) + 2.8. This expression is almost independent of the 
number of branons. From the numerical values of Tj in terms of /, it is pos- 
sible to obtain Qsrh 2 from (19). In this case we have considered two kinds of 
limits. On one hand those coming from the total dark matter of the universe 
flBrh 2 < 0.129 — 0.095 in Fig. 2. On the other hand, more constraining lim- 
its on the hot dark matter energy density can be derived from a combined 
analysis of the data from WMAP, CBI, ACBAR, 2dF and Lyman a [1]. The 
bound reads Q Br h 2 < 0.0076 at the 95%C.L. and it is obtained thanks to the 
fact that hot dark matter is able to cluster on large scales but free-streaming 
reduces the power on small scales, changing the shape of the matter power 
spectrum. In Fig. 3 we have plotted the corresponding limits in the f — M 
plane. Notice the abrupt jump around / ~ 60 GeV in Figs. 2 and 3 (and also 
in Fig. 4). This / value corresponds to a decoupling temperature of T ~ 170 
MeV which is the assumed value for the QCD phase transition. Thus, this 
jump is due to the sudden growth in the number of effective degrees of free- 
dom when passing from the hadronic to the quark-gluon plasma phase. The 
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exclusion areas depend on the number of branon species and we have plot- 
ted them for AT = 1, 2, 3, 7 in Figs. 2,3. The validity of the previous limits 
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Figure 2: Exclusion plot for hot relics in models with N = 1, 2, 3, 7 branons. 
For a given N, the shaded area corresponds to the total dark matter limit 
flBrh 2 = 0.129 — 0.095, therefore the region above such area is excluded by 
branon overproduction. 

requires that branons were relativistic particles at freeze-out. Therefore we 
require < 3, which implies that the bounds do not work for / < 10 -4 
GeV. The curve Xf = 3 in the hot relic case is also plotted in Fig. 6. 

As commented in the introduction, in all the previous calculations we are 
assuming, apart from / <C Mp, that the evolution of the universe is standard 
up to a temperature around /. In fact, the effective Lagrangian (4) is only 
valid at low energies relative to / and therefore, it is this scale what fixes 
the range of validity of the results. We have checked that our calculations 
are consistent with these assumptions since the decoupling temperatures are 
always smaller than / in the allowed regions in Fig. 6. 

6 Bounds from nucleosynthesis 

As commented above, one of the most successful predictions of the standard 
cosmological model is the relative abundance of the light elements. The cal- 
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Figure 3: Exclusion plot for hot relics in models with N = 1, 2, 3, 7 branons. 
For a given N, the curve corresponds to the hot dark matter limit flBrh 2 = 
0.0076, therefore the region above such curve is excluded by hot branon 
overproduction. 



culated abundances are very sensitive to certain cosmological parameters, in 
particular it has been shown that the production of 4 He increases with in- 
creasing rate of expansion H. From (18) we see that the Hubble parameter 
depends on the effective number of relativistic degrees of freedom g e ff- Usu- 
ally, this number is parametrized in terms of the effective number of neutrino 
species N u = 3 + AN U as: 

g eff (T ~ MeV) = g s ef M f + gffi < 10.75 + -^N u (34) 

where T ~ MeV corresponds to the universe temperature during nucleosyn- 
thesis. In the SM with three massless neutrino families, we have gf^j{T ~ 
MeV) = 10.75 corresponding to the photon field, the three neutrinos and 
the electron field. In (34), in order to avoid deviations of the predicted 
abundances from observations, the conservative limit AN U = 1 for the con- 
tribution from new physics is usually imposed [17], i.e. there could be only 
one new type of light neutrino. 

Including branons, the number of relativistic degrees of freedom at a given 
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temperature T is given by: 



i 



9eff(T)=g™(T) + N^) (35) 

where gf^{T) is the contribution from the SM particles, T B denotes the tem- 
perature of the cosmic branon brackground and we are assuming that there 
are no additional new particles. If branons are not decoupled at a temper- 
ature T then Tb = T, i.e. they have the same temperature as the photons. 
On the other hand, if they are already decoupled then its temperature will 
be in general lower than that of the photons. In order to calculate it, we use 
the fact that the universe expansion is adiabatic. Let us write: 

h eff (T) = h s e f f (T) + N (^) 3 (36) 

where hffj(T) includes only the contributions from SM particles. If at some 
time between branon freeze-out and nucleosynthesis, some other particle 
species becomes non-relativistic while still in thermal equilibrium with the 
photon background, then its entropy is transferred to the photons, but not 
to the branons which are already decoupled. Thus, the entropy transfer in- 
creases the photon temperature relative to the branon temperature. The 
total entropy of particles in equilibrium with the photons remains constant 
i.e.: 

h e e q ff a 3 T 3 = constant (37) 

and since the number of relativistic degrees of freedom h e ^j has decreased, 
then T should increase with respect to T B . Thus, we find: 

where Tf >B is the branon freeze-out temperature and we have used the fact 
that for particles in equilibrium with the photons g^jj = h e e q ^. Using (35), 
we can set the following limit on the number of massless branon species iV: 

-AN >n( -^-V = n( ^ /(T "" c) (39} 

If branons decouple after nucleosynthesis we get the direct limit: 

N < ^AN U (40) 
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If they decouple before, we have g e ^^{T n 
ingly, we can rewrite the bound as: 



10.75, as seen before. Accord- 



N < 7 -AN v 



'9eff(Tf, 



B 



4/3 



10.75 



(41) 



Taking AN U = 1, the relation between the freeze-out temperature Tj s and 
the brane tension scale / that we have obtained for hot relics can be used to 
get limits on the number of branons N (see Fig. 4). Thus, for / < 10 GeV, 
we get N < 1. This result is obtained using eq.(40) for / < 3 GeV (which 
corresponds to T^b ^ 1 MeV) or eq.(41) otherwise. However the limits are 
less restrictive in the range / ~ 10 — 60 GeV. In this case we get N < 3. 
Above the QCD phase transition which, as commented before, corresponds 
to / ~ 60 GeV, the bound rises so much that the restrictions are very weak. 
Notice that we are taking gfHiT > 300 GeV) = 106.75, i.e. we only include 
the minimal Standard Model matter content with a Higgs doublet and three 
massless neutrinos. 
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Figure 4: Restrictions from nucleosynthesis on the number of massless branon 
species N SIS db function of / for AN V = 1. 

Concerning the value of AN U , more constraining analysis using only BBN 
suggest AN U = 0.5 [19]. However, using also WMAP results, the constraints 
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are N v = 2. (95% C.L.), which are only marginally consistent with 
the LEP measurements of the number of neutrino families [2]. Such values 
would severely constrain the number of any new light particles present during 
nucleosynthesis. However, it has been suggested that these results could 
be potentially affected by systematics errors in the BBN predictions of the 
primordial abundances [2, 20]. 

7 Bounds from supernova SN1987A 

Important astrophysical bounds on the brane tension scale can be obtained 
from the energy loss in supernovae [12]. Such energy loss is carried away 
essentially by light particles, i.e. photons and neutrinos if we restrict our- 
selves to SM particles only. However, if the branon mass is low enough, we 
expect branons to carry some fraction of the energy, whose importance will 
depend on their couplings to the SM particles. In this section, we study the 
constraints on / and M imposed by the cooling process of the neutron star 
in supernovae explosions. We will perform an analysis of the energy emission 
rate from the supernova core, similar to that done by Kugo and Yoshioka 
for massless branons [12]. The aim is to extend their study to arbitrary 
M and compare with the colliders and cosmological bounds. In their work, 
they consider the channel corresponding to electron-positron pair annihila- 
tion. Although this contribution could be subdominant, it will allow us to 
get an order of magnitude estimation of the branon effect. 

If branons are produced in the core, they can be scattered or absorbed 
again depending on their couplings to the SM particles. Only if the bra- 
non mean free path L inside the neutron star is larger than the star size 
(R ~ O(10) Km), they could escape and carry the energy away. For a mas- 
sive branon (M 3> T$n), we have L ~ (87r/ 8 )/(M 2 T| 7V n e ), where n e is the 
electron number density in the star. Therefore the restrictions we will obtain 
will be valid typically only for / > 5 GeV. These restrictions appear due 
to the fact that the emitted energy in the form of branons could spoil the 
agreement between the predictions for the neutrino fluxes from supernova 
1987A and the observations in the Kamiokande II [21] and 1MB [22] detec- 
tors. Branons could shorten the duration of the neutrino signal if the energy 
loss rate per unit time and volume is Q > 5 x 1CT 30 GeV 5 . In particular, 
the contribution of the mentioned channel to the volume emissivity has the 
form: 
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Figure 5: Exclusion regions from supernovae cooling by branons for T$n — 
10, 50, 70 MeV and N — 1. The solid lines come from the limits on the 
volume emissivity, and the dashed lines are the L — 10 Km limits on the 
branon mean free path. 



r 2 ( d 3 k- 1 
Q Br (f,M) = J I] I {27T) 32 El 2f j {El + E2)2g ^+e--^(a, /, AO (42) 

where i refers to the electron (1) and positron (2) particles, whose masses 
can be neglected in the supernova core. The chemical potential in the Fermi- 
Dirac distribution function /, = 1/ ^ Ei ' T ~ ti ' T ^ + 1) can be estimated as 
/i ~ (37i 2 n e ) 1 ^ 3 with the number density of electrons: n e ~ 1.4 x 10~ 3 GeV 3 . 
With these assumptions, Qb t is given by: 



poo poo rl— ZM 

Qbt = / dE x dE 2 d(cos){E 1 + E 2/ 

Jo Jm2/e 1 J -i 

N[2E 1 E 2 (2E 1 E 2 (1 - cos) - AM 2 )] 5 / 2 (1 - cos) 3/2 



(2 7 r) 5 7680r(i + e £ ^) (l + e^) ^ 

In fact, it is possible to calculate analytically the angular integral, whereas 
the integral over the two energies has been performed numerically. The 
corresponding constraints depend on the supernova temperature (T$n) and 



19 



the number of branons (N). In Fig. 5 we show the limits on / and M for 
T SN = 30, 50, 70 MeV and N = 1. 

It is interesting to note that for a branon mass of the order of the GeV, 
the restrictions on the brane tension scale disappear even for T S n = 70 MeV 
due to the limitations in the mean free path discussed above. 

8 Conclusions 

In this work we have studied the limits that cosmology and astrophysics 
impose on the brane-world scenario through the effects of massive branons. 
Using the effective low-energy Lagrangian for massive branons interacting 
with the SM fields, we have computed the annihilation cross sections of bra- 
non pairs into SM particles. From the solutions of the Boltzmann equation 
in an expanding universe, we have studied the freeze-out mechanism for bra- 
nons and obtained their thermal relic abundances both for the cold and hot 
cases. Comparing the results with the recent observational limits on the total 
and hot dark matter energy densities, we have obtained exclusion plots in the 
f — M plane. Such plots are compared with the limits coming from collider 
experiments and show that there are essentially two allowed regions in Fig. 
6: one with low branon masses and large brane tensions (weak couplings) 
corresponding to hot branons, and a second region with large masses and 
low tensions (strong couplings) in which branons behave as cold relics. In 
addition, there is an intermediate region where / is comparable to M, which 
is precisely the region studied in [13], and where branons could account for 
the measured cosmological dark matter. 

Using the nucleosynthesis limits on the number of relativistic species, we 
set a bound on the number of light branons in terms of the brane tension. We 
see that if branons decouple after the QCD phase transition corresponding to 
/ < 60 GeV, the limits can be rather stringent N < 3, whereas they become 
very weak otherwise. 

Finally, we have analysed the possibility that massive branons could con- 
tribute to the cooling of a supernova core. After estimating the energy loss 
rate, we again get some limits on the / and M parameters which are com- 
pared to the previous ones. It is shown that they are not competitive with 
those coming from LEP-II. 

In conclusion, cosmology imposes limits on the BWS which are comple- 
mentary to those coming from collider experiments and astrophysics. Al- 
though the combination of both bounds excludes an important region of the 
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Figure 6: Combined exclusion regions in a model with a single branon from 
total and hot dark matter, LEP-II single photon events [11] and supernovae 
cooling. The (blue) solid line on the right correspond to the cold dark matter 
limit. The two (red) solid lines on the left correspond to hot dark matter: the 
thicker one comes from the total dark matter range VLs r h 2 = 0.129 — 0.095, 
whereas the thin one is the hot dark matter limit Qsrh 2 = 0.0076. The two 
dashed lines correspond to Xf = 3 for hot (upper line) and cold (lower line) 
dark matter. 



parameters space, still there are brane-world model which could be compati- 
ble with observations. Future hadronic colliders such as LHC, Tevatron-II or 
the planned linear electron-positron colliders, and the possibility of detect- 
ing dark matter branons directly [13] or indirectly will allow us to explore a 
wider region of the parameter space. Work is in progress in these directions. 
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A Branon production and annihilation cross 
sections with SM particles 

In this section we present the branon production and annihilation cross sec- 
tions in processes involving SM particles. The results are presented with all 
the internal degrees of freedom summed for the final particles and averaged 
for the initial ones. We have used the Feynman rules given in [11], where N 
is the number of branons. 

A.l Scalar s 

<7i: $t(pi),$(p 2 ) — >7r(P3),7r(p 4 ) 



N 



o x = 



4M2 )[_ s ( 8m | + s )( s _4M 2 ) 2 



7680/%s^ (s-4m|)' 
+ (2m| + s) 2 (23M 4 - UM 2 s + 3s 2 )}. (44) 

a 2 : 7r(pi),7r(p 2 ) — > ^(ps),^) 



<?2 



' s 4m *)[_ s ( 8m | + S )( S _4M 2 ) 2 



3840N f 8 ns\ (s - 4M 2 ) 
+ (2m| + s) 2 (23M 4 - 14M 2 s + 3s 2 )]. (45) 

The scalar results are given for a complex scalar such as charged pions 
or kaons. For a real scalar like the Higgs field or the neutral pion, the 
only change comes from the branon annihilation cross section that should be 
divided by two. 

A. 2 Fermions 

^3- *P + (Pl),1p~(P2) ► 7r(P3),7T(p 4 ) 

-3 = v/( S - 4M 2 )( S - 4m 2 )[(s - 4M 2 ) 2 

2m 2 

+ ^(23M 4 - 14M 2 s + 3s 2 )]. (46) 
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a 4 : 7r(pi),7r(p 2 ) — >ip + (p 3 ),ip (p 4 ) 



(7 4 = 



:S - 4m ' )3 [( S -4M 2 ) 2 



3840^/%^ (s-4M 2 ) 
2rr? 2 

+ ^(23M 4 -14M 2 s + 3s 2 )]. (47) 



These results are valid for a Dirac fermion of mass m^. For massless 
Weyl fermions, we should multiply the branon production cross section by 
two, whereas the annihilation cross section should be divided by the same 
factor of two. 

A. 3 Photons 

<^5- 7(Pl)»7(P2) ► 7r(P3),7r(P4) 



(7 6 : 7r(pi),7r(p 2 ) — >7(Ps),7(P4) 



1 s(s-4M 2 ) 2 
° % ~ 1920iV/ 8 7r ^_ 4M£ • (49) 



A.4 Z 



cr 7 : Z{j>x),Z{pi) — ► 7r(p 3 ),7r(p4) 



TV /s-4M 2 9 9x9 

= 69lW%JV^iM| |3S<8M| + S)(S - 4M) 

+(12M§ + 4sM| + s 2 )(23M 4 - 14M 2 s + 3s 2 )]. (50) 

(7 8 : 7r(pi),7r(p 2 ) — ► Z(p 3 ),Z(p A ) 



+(12M| + 4sM| + s 2 )(23M 4 - 14M 2 s + 3s 2 )]. (51) 
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A.5 W ± 



a 9 : W ± (pi),W^(p 2 ) ^ 7r(p 3 ),7i(p 4 ) 



N s- AM 2 r . 9 . . 9x9 

ff9 = emw^V^«¥ [3s(8M " + s)(s - 4M) 

+ (12M^ + 4sM^ + s 2 )(23M 4 - UM 2 s + 3s 2 )]. (52) 
a w : 7r(pi),7r(p 2 ) — ► W^Cps), W T (P4) 



+ (12M^ + AsM^ + s 2 )(23M 4 - 14M 2 s + 3s 2 )]. (53) 



A. 6 Gluons 

<^n- g(Pi),g(P2) — ► 7r(p 3 ),7r(p 4 ) 



N AM 2 9n9 

ff » = eiSoT^V 1 -— S(S " 4M) ' (54) 
0-12: 7r(pi),7r(p2) — > g(p3),g(p*) 



C12 = 



1 ^-4M 2 ) 2 
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